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Abstract. 

An explicit formula for the resistance between two nodes in a network with a non- 
symmetric Laplacian matrix L is obtained. This is of great advantage e.g. in 
electronic circuit fault analysis, where non- linear systems have to be solved repeatedly. 
Analysis time can be greatly reduced by utilization of the obtained formula. The 
presented approach is based on the "mutual orthogonality" of the full system of left 
and right-hand eigenvectors of a diagonalizable matrix L. Some examples are given 
that demonstrate the accuracy of the approach on circuit networks. 
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1. Introduction 

The computation of the resistance between two nodes in a resistor network is a classical 
problem in electric circuit theory, indeed it was studied by numerous authors over many 
years. Since electrical potentials on a grid are governed by the same difference equations 
as those occurring in other problems, the computation of the resistances can be applied 
to a wide range of problems, such as lattice Green's functions, random walks, the theory 
of harmonic functions and others. Due to variety of applications of the resistance 
problem and various methods of solution, this problem is still actual and intensively 
studied [T]. 

The resistance problem was often transformed in that of solving the difference 
equations, which is most conveniently carried out for infinite networks. Little attention 
was paid to finite networks, though these are occurring in real life. In 2004, an important 
result for finite lattices was published in [2]. He expressed the resistance between two 
arbitrary nodes in a resistor network in terms of eigenvalues and eigenvectors of the 
Laplacian matrix L associated with the network. In [2], the Laplacian L is a real 
symmetric N x N matrix with N real eigenvalues Aj and an orthonormal basis of 
eigenvectors i/'i = (^a, • • • , ipiN)* ■ The resistance between the nodes a and (3 is given 
by 

i=2 1 

However, if the Laplacian matrix of the network is non-symmetric, the relation ([T]) 

fails. 

In this work, we study resistance networks with a real Laplacian matrix L, 
symmetric or not. The Laplacian L is a real diagonalizable N x N matrix with N real 
eigenvalues Aj and two basis of "mutually orthogonal" left and right-hand eigenvectors 
if* = ((fa, . . . , ipiN) and if>i = (tpa, . . . , To our best knowledge, this is the first 

work that presents an explicit formula for the total resistance between two nodes in a 
resistor network with a real Laplacian matrix, in general non-symmetric. The obtained 
formula 

N 1 

is an extension of Q. The method used in this paper is similar to that of [2J, where 
the orthonormal basis of eigenvectors of the matrix L is replaced by two "mutually 
orthogonal" basis of left and right-hand eigenvectors of L. 

The organization of the paper is as follows. In section [2] our motivation coming 
from the fault analysis of integrated circuits is presented and the nodal conductance 
matrix of a circuit, which is the Laplacian L, and some of its properties are described. 
Section [3] gives a theoretical framework, which is based on the theory of eigenvalues and 
eigenvectors of real matrices. The existence of "mutually orthogonal" basis . . . , ip* N 
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and ipi, . . . , ipN of left and right-hand eigenvectors of a real diagonalizable matrix A is 
proved. Keeping the vectors ip\, . . . , ip* N intact, a convenient basis of right eigenvectors 
is constructed. Then, the main result is formulated as theorem and proved in the second 
part of section |3j Section [4] is devoted to calculating the two-point resistance for several 
circuit networks. The formula proved in section [3] is demonstrated on a simple bridge for 
two cases: simple eigenvalues and multiple eigenvalues. Also, an example of resistance 
calculation for a real life circuit, an operational amplifier, is given. 

2. Fault analysis - motivation 

Today's integrated circuits (ICs) represent complex systems realized on a small area 
using advanced fabrication technologies. At the early age of monolithic integration (early 
sixties), the fabrication process was dominated by yield as low as 10% due to poor process 
control [3]. As the fabrication process has evolved, the process control reached acceptable 
limits and the yield improved. On the other hand, with the increasing resolution, the 
fabrication process became more sensitive to impurities of the environment and had to 
be moved to controlled areas, where the ambient impurity level was under control. 

These days, a typical mixed-signal IC can consist of thousands to millions of 
transistors, occupy an area of several tens of square mm, while the fabrication process 
yield with resolution as high as tens of nanometers is well above 90%. Anyhow, even in 
a well controlled environment, the ambient still exhibits certain level of impurities. 
Besides, the fabrication process control is becoming even more complex and thus, 
resulting in a necessity to test the produced ICs. 

During the test development process, so called fault simulations have to be 
conducted to determine the fault coverage of a given test. If performed in a traditional 
way, this step might be rather time consuming, since it requires a huge amount of 
simulations to be executed due to large amount of circuit nodes and possible faults 
taken into account. During the simulation, a set of nonlinear equations describing the 
circuit has to be solved using an iterative algorithm, e.g. Newton-Raphson. Application 
of such an algorithm to circuit analysis is described in [8]. 

The numerical approach to DC fault simulations, presented in this paper, evaluates 
the resistance between two arbitrary nodes of the circuit in a more time efficient way. 
This evaluation requires only one simulation for the given circuit configuration, providing 
information on the operating point, to be executed. It means that a system of nonlinear 
equations describing the circuit has to be solved only once for the given operating point. 
From this operating point information, the circuit's nodal conductance matrix can be 
created and the resistance between the selected two nodes can be calculated. As has 
been already shown for symmetrical matrices, such an approach can lead to a significant 
analysis time reduction [3] . This resistance calculation can be used in the development of 
resistance-based tests, e.g. Iddq testing, where the resistance is indirectly evaluated by 
measuring the current flowing from the supply source [51 E]- For a circuit network, where 
active elements such as MOS transistors are present, the nodal conductance matrix can 
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become non-symmetrical, since voltage-controlled current sources (VCCS) are used to 
model such devices, as described below. 

2.1. Nodal conductance matrix 

Formulation of the nodal conductance matrix, which is the Laplacian, is based on the 
Kirchhoff current law (KCL) which states the following [7j 

N 

J2' c n(yi-Vj) = Ii t = l,2,...,N. (3) 

i=i 

Here, is the conductance of a resistor connected to nodes i and j, Vi and Vj are the 
electrical potentials at nodes i and j, and is the current flowing into the network 
at node i. The prime denotes the omission of the term i = j. As a consequence, the 
following constraint applies: 

N 

£ J « = °- (4) 

1=1 

In general, a circuit of iV nodes can be described by a matrix equation of the form: 

LV = f (5) 

or 



/ Cn C12 ... C 1N \ /Vi\ ( h \ 
C21 C 2 2 ■ ■ ■ C 2 N V 2 

\cni c N2 ... c NN J \V N J \l N J 



(6) 



where L is the Laplacian (nodal conductance matrix) with elements cn, . . . , cnn, V is 
a column vector of node potentials, and / is a column vector of independent current 
sources [S]. 

To compute the resistance between two arbitrary nodes a and /3 of the network, an 
external source has to be connected to these nodes. Let the current flowing from the 
source be / and the potentials at the two nodes V a and Vp. Then the desired resistance 
is determined by the Ohm's law 



R n 



V n -V, 



(7) 



Now, the computation is reduced to solving (|5j for V a and Vp with the current given by 
I i = I{5 ia -5 i p). (8) 
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An important property of the L matrix is 



A' 



N 



(9) 

i=l i=l 

From this equation follows, that since the KCL requires the sum on the right side to be 
zero, the sum of conductances in each column of L is equal to zero. Another important 



property is given by (10): 

N 



0. 



(10) 



Since Vj is a nonzero voltage source, the sum must be equal to zero. Therefore, the sum 
of conductances in each row of L is equal to zero. 

As long as the network consists of resistances only, the Laplacian remains 
symmetric, since for a resistive element of a conductance c connected between the nodes 
i and j, the nodal conductance matrix L R can be expressed as follows: 




If this matrix is added to the Laplacian of the network, the symmetry of the Laplacian is 
maintained. On the other hand, if a VCCS is present in the network, the Laplacian can 
become non-symmetrical in cases when the input terminals are connected to different 
nodes as the output ones. In general, the nodal conductance matrix Lyccs of a VCCS, 
as depicted in figure [TJ can be written as 



'VCCS 
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(12) 




Figure 1: Voltage controlled current source 

Since most of the real life circuits contain active elements that are represented by 
VCCSs, an extension of the approach published in [2| to non-symmetrical Laplacians is 
of great importance. 
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3. The two-point resistance theorem 

In various applications based on matrix theory, a diagonalization of a square matrix is 
needed. A real N x N matrix A is diagonalizable, if there exists an invertible matrix P 
such that P X AP = D, where D is a diagonal matrix. Here, column vectors of P are 
right-hand eigenvectors of A and row vectors of P _1 are left-hand eigenvectors. 

When the matrix A is symmetric, the eigenvectors can be chosen to form an 
orthonormal basis of M. . In this case, P is an orthogonal matrix and P _1 equals 
the transpose of P. The diagonal elements of D are eigenvalues of both A as D. 

When a diagonalizable matrix A is non-symmetric, the orthogonality of eigenvectors 
fails. Nevertheless, since neither the orthogonality of right eigenvectors nor that of left 
ones is necessary in our approach, this problem can be overpassed. What we need, is 
the "mutual orthogonality" of the left and right eigenvectors. Then, if $ and \I/ denote 
matrices of left and right-hand eigenvectors of a square matrix A, the matrix $A^ is 
diagonal. Its diagonal elements are not necessarilly the eigenvalues of A. 

We state our main result as a theorem. In its proof, we use a similar method as in 
the proof of Theorem in [2J . 

From Kirchhoff's law we know that at least one eigenvalue of Laplacian matrix L 
is zero. We consider this eigenvalue simple and put Ai = 0. The associated left and 
right-hand eigenvectors are <pi = (1, 1, . . . , 1) = The condition of orthogonality of 
left and right-hand eigenvectors (<p*ipj = 0) (i ^ j) can be fulfilled for arbitrary 
diagonalizable matrix that satisfies any other condition of the following theorem. 

Theorem. Consider a resistance network with N nodes defined by a diagonalizable 
Laplacian matrix L not necessarily symmetric. Suppose that L has N real eigenvalues 
Aj (including the multiplicity) with left and right-hand eigenvectors ip* = ((fin, . . . ,<Pin) 
and if>i = (ipn, . . . ,ipm)*, where (fi and iftj are orthogonal if and only if i ^ j. Denote 
Ai the zero eigenvalue. Then the resistance between nodes a and (5 is given by 

N 1 

Proof. Since one of eigenvalues of L is zero, introducing of eL to the Laplacian L will 
be helpful. Thus we will manipulate the Laplacian L(e) = L + eL. The eigenvalues of 
L(e) are e, A2 + e, . . . , Aat + e and for a suitable value of e, no one of them equals zero. 
The left and right-hand eigenvectors are those of L. Similarly as in [2], we will compute 
the inverse function G(e) of L(e). 

Naturally, L(e) is not symmetric, because L is not. Non-symmetric matrices L and 
L(e) are diagonalizable as follows: 

$L* = A and $L(e)* = A(e), (13) 

where $ is the matrix of left-hand eigenvectors (rows) and ^ that of right-hand 
eigenvectors (columns) of L and L(e). Finally, A(e) is a diagonal matrix whose diagonal 
elements are (Aj + e)(p*if>j. 
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Both and ^ are invertible, because {ipi, . . . , 95 at} as well as {if>i, . . . , i^n} define 
a basis of M. N . From (13) we obtain 



*- A L- A (e)*- A = A" A (e) 



Hence 



G(e) = *A _1 (e)*. 

Really, G(e) is an inverse matrix of L(e), because 

G(e)L(e) = (*A- 1 (e)$)($- 1 A(e)*- 1 ) = I, 

as well as 

L(e)G(e) = ($- 1 A(e)*- 1 )(*A- 1 (e)$) = I. 



(14) 



(15) 



From (15) it follows that elements of the matrix G(e) are 

G "W = E*.( (A| + l )V fJ^ = ^+E 



i/3 



(16) 



So, the function G(e) is explicitly expressed. 

Now, return to the equation L(e)V(e) = I and multiply it from the left by G(e) 
We obtain V(e) = G(e)J. Coordinates of the vector V(e) are 



N 



Vi{e) = Y,Gij(e)l3 for i = l,2,...,JV. 

i=l 



(17) 



Denote in (16) 



A? 



M = E 



i=2 



(\i + e)cp*il>i' 



The substitution of (e) = (iVe) 1 + (e) in QTtT) yields 



N ( 1 \ 1 * " 

= E ( jvg + ) J i = jvg E J i + E 9aWr 
3=1 v J 3=1 3=1 



(19) 



A? 



A 



Since 7j = 0, we have V*(e) = gij(e)Ij and we can put e = O.This gives 



i=l 



i=i 



A 



^ = E^> 



(20) 



i=i 



where V* and gry replace V^(0) and <7y(0). 
We use g in = I -1 ^ - V» which is @ 



1 / N N \ N L 

J E 9°* I i - E I i = E ~ 9Pj)j- 
\i=i i=i / j=i 



(21) 
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Due to (1st) we obtain 



N 

Ra/3 = 2J (g a j — 9/3j)($ja — 5jp)- 
3=1 

The last expression can be simplified as 

Ra/3 = (gaa ~ 9/3a)(l ~ 0) + (#q,/3 - 5 , /3/?)(0 ~ 1). 



We introduce (18) with e = and (23) becomes 

AT N 



R, 



aj3 



E 



N N 



E 



A simple rearrangement of the last relation returns (|2]). 



(22) 
(23) 

(24) 

□ 



4. Examples 

In the following examples, the proposed resistance calculation approach is demonstrated 
on simple networks (circuits). The calculations were conducted using GNU Octave lan- 
guage. The circuit analysis results were obtained from Cadence IC design environment. 



Example 1 (simple eigenvalues). 



Figure 2: A network of 4 nodes with VCCS 

Consider the network depicted in figure [2] For this network, the nodal conductance 
matrix, which is the Laplacian L, is the following non-symmetric matrix, where 
C\ = 1/Ri and c 2 = 1/-R 2 : 



/ 2ci — ci — cA 

— ci + G 2ci + c 2 — ci — c 2 - G 

— ci 2ci — ci 

\— ci - G — c 2 — ci 2ci + c 2 + G/ 



The eigenvalues are Ai = 0,A 2 = 2ci,A3 = 4ci,A4 = 2ci + 2c 2 + G. Since 
the eigenvectors are not required to have the norm 1, we take those with the least 
complicated entries. The left and right-hand eigenvectors are 
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^2 


^3 




( 1 


1 


1 


G(a - 2c 2 - G)\ 


1 





-1 


( Cl -c 2 )(2c 2 + G) 


1 


-1 


1 


-dG 


V i 





-1 (2c 2 - 


\-G)(- Cl + c 2 + G)J 





ll>2 




^3 




V>4 


( 1 


-2c 2 - G 


-2 Cl 


+ 2c 2 


+ G 




1 


G 


2Ci- 


-2c 2 - 


-3G 


1 


1 


2c 2 + G 


-2ci 


+ 2c 2 


+ G 





\ 1 


-G 


2ci 


-2c 2 


+ G 


-1/ 



Using (J2J), the resistance between two arbitrary nodes of the network can now be 
calculated. Below, examples of resistance calculation follow. 

1 1 



R 



13 



(y?31 - ^33) ("031 - ^33) + 



R 



1 



(V?41 - V43) (^41 - V> 



43 ^ 



21 



\ 2 <P* 2 ll) 2 
1 

+ 



((/? 2 2 - ^24)(^22 - ^24) + 
-(<^42 - <^44) (^42 - ^44) 



1 



"(^32 - ^34) (^32 - ^34) + 

1R\R 2 



A 4 ^4 " y 2R 1 + 2R 2 + GR 1 R 2 

Set the values of parameters i?i = 200f2, i? 2 = 4000fi, G = 0.035. All the possible two- 
node resistances constitute aniVxiV upper triangular matrix where N is the network 
node count. The matrices of calculated resistance R c and resistance from circuit analysis 



R a are shown below. One can observe that R c 



Ra- 



Rr 



/0 1.864E + 02 2.000E + 02 3.827E + 01\ 

1.124E7 + 02 4.938E + 01 

1.124E + 02 

\0 0/ 



1.864E + 02 2.000E + 02 3.827E + 01\ 

1.124E + 02 4.938£ + 01 

1.124E + 02 

0/ 

Example 2 (multiple eigenvalues). 

Consider the network depicted in figure |3} Let c\ = 1/R. For this network, the 
following non-symmetrical nodal conductance matrix, which is the Laplacian L, can be 
created: 



R„ 






Vo 
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Figure 3: A network of 4 nodes with VCCS 



/ 3ci -ci -ci -cA 

— Ci 3ci — ci — Ci 

-ci + G — ci 3ci —c\ — G 

\—c\ — G — c\ — ci 3ci + Gy 



The eigenvalues are Ai = 0, A 2 = 4ci,A 3 = 4ci,A 4 = 4ci + G. Note the multiple 
eigenvalue A2 = A3 = 4ci. Left and right-hand eigenvectors <pi and ^ of this matrix 
are as follows: 



ifil y> 2 Vs ip A 



( 1 
1 

1 

V 1 



1 

-1 









-1/ 



( 1 
1 

1 

V 1 



^2 */>3 
1 1 



1 

-3 
1 



i>4 

°\ 



1 

-1/ 



As can be observed, the left and right-hand eigenvectors are independent of C\ and G. 
An example calculation of the first row of the resistance matrix R c is given below. We 
can use substitution c\ = 1/R for this calculation. 

1 R 1 R 1 R 

A 2 <^2?/>2 16' 



A 3 ^3 8' A 4 ^4 4 + Gi?' 



R12 — 



R 
16' 



1 + 1)(1-1) + |(1 + 1)(1 + 1) + 



R 



A + GR 



l-0)(0-0) 



+ -4 + = -, 
8 2 ' 



*is = - ^(-1 - 1)(1 + 3) + f (1 - 0)(1 + 1) + 4^(1 " 0)(0 - 1) 



R R R 
2 + J~ 4 + GR 



p/3 1 
!, 4 4 + Gi? 



^14 = - 1)(1 - 1) + f (1 - 0)(1 - 1) + + !)(° + !) 



+ + 



2R 
4 + GR' 
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Now, set the values of parameters R = lfi and G = 2S. The values of resistances are 
R12 = 1/2, #13 = 7/12 and i? i4 = 1/3. The complete resulting matrices of calculated 
resistance R c and resistance from circuit analysis R a are shown below. One can observe 
that R c = R a . 



R f 





1/2 


7/12 


l/3\ 




/o 


5.000E - 01 


5.833£ - 01 


3.333E 


-01\ 








1/2 


5/12 


Ra — 








5.000E - 01 


4.167E 


-01 











1/3 













3.333£ 


-01 


\o 












\o 










0/ 



Example 3 (operational amplifier). 

The circuit depicted in figure [4] represents a typical topology of an operational 
amplifier that is one of the most frequently used analog building block. 



ah 



o 



4T 



gnd 



4 



n2 



rrr\ 

UMl M2 r 1 



3h 



T 



1 



n4 



Figure 4: Circuit schematic of an operational amplifier 



This particular circuit has been designed for a standard 0.35/im mixed-signal CMOS 
process. Parameters of the devices obtained from a typical .op (operating point) analysis 
are shown in table Q] 

Table 1: Device parameters of the operational amplifier 



gm[S] gds[S] Ron[ft] 

Ml 5.155032^- 04 2.816268S - 03 - 

M2 5.155032£- 04 2.816268S - 03 - 

M3 5.953932£ - 04 3.403570S - 06 - 

M4 5.953932£ - 04 3.403570S - 06 - 

M5 4.195651£- 03 9.344114S - 06 - 

M6 5.247305£ - 04 2.065995S - 06 - 

Rl - - 6.847545^ + 04 



These parameters are used to create an electrical model of the circuit, where only 
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resistances and voltage-controlled current sources (VCCS) are used. MOS transistor is 
modeled by a VCCS, which is represented by a transconductance gm and a resistor of 
value 1/gds connected between its output terminals, as depicted in figure [51 



D 



1/gds 



s 



Figure 5: Simplified electrical model of a MOS transistor 

Using the parameters from table [TJ the following non-symmetrical nodal admittance 
matrix L of this circuit can be created for the internal nodes: 



5AUE 


-04 


0.000E 


+ 00 


0.000E + 00 


0.000E + 00 


-5AUE 


-04\ 


4.196E 


-03 


6.673E 


-03 


-2.816E 


-03 


-2.816£-03 


-5.236E 


-03 


0.000E 


+ 00 


-3.332£ 


-03 


3.415E 


-03 


0.000E + 00 


-8.329E 


-05 


0.000E + 00 


-3.332£ 


-03 


5.954E 


-04 


2.820£ - 03 


-8.329£ 


-05 


4.737E 


-03 


-9.344E 


-06 


-1.194E 


-03 


-3.404E - 06 


5.9UE 


-03/ 



The node Vdd is excluded from this matrix, since when assuming the impedance 
between two nodes, all the voltage and current sources have to be substituted by their 
equivalent circuits. Since no parasitic resistance is taken into account, a voltage source 
is substituted by a short, while a current source is represented by an open. In our case, 
it means that the node Vdd becomes shorted to node gnd. The order of network nodes 
assumed for this matrix is nl, n2, n3, n4, gnd. 

For the nodal admittance matrix, which is the Laplacian L, the eigenvalues obtained 
are A a = 0, A 2 = 1.617872E - 04, A 3 = 2.819672E - 03, A 4 = 8.008586E - 03, A 5 = 
8.402968£ - 03. 

The corresponding left and right-hand eigenvectors (in columns) of the matrix L 
are as follows: 
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/0.447214 
0.447214 
0.447214 
0.447214 

V0.447214 



¥>2 

0.890916 
-0.250715 
-0.223688 
-0.265848 
-0.150665 



vl 




-0.707107 
0.707107 




0.606347 
0.233943 
0.034793 
-0.126481 
-0.748603 



^5 

-0.596298^ 
-0.286327 
0.002639 
0.143977 
0.736009/ 



/0.447214 
0.447214 
0.447214 
0.447214 

VO. 4472 14 



V'2 

-0.638056 
0.363601 
0.360920 
0.360920 

-0.447385 



^3 

-0.038249 
0.098610 
0.574328 

-0.795646 
0.160957 



tj>A 

0.023987 
-0.654286 
0.480565 
0.480565 
-0.330831 



V>5 

-0.017999^ 

0.698457 
-0.470914 
-0.470914 
0.261369/ 



The resulting matrices of the numerically calculated resistance R c and the resistance 
from the circuit analysis R a are shown below. One can observe that R c = R a . 

/0 9.276302E + 03 9.098426E + 03 9.448754E + 03 1.847062E + 03\ 

3.000575E + 02 3.546762E + 02 8.490922E + 02 
= n n n 7.093024E + 02 



Vo 









8.317060E + 02 
1.182034E7 + 03 
0/ 







Vo 



9.276302E + 03 







9.098426E + 03 
3.000575E + 02 






9.448754E + 03 
3.546762E + 02 
7.093024E + 02 





1.847062E + 03\ 
8.490922E + 02 
8.317060E + 02 
1.182034^ + 03 
0/ 
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5. Conclusion 

A formula for resistance computation between two arbitrary nodes of a network with 
non-symmetric Laplacian matrix L has been formulated and proved in this paper. The 
calculation is based on eigenvalues and eigenvectors of the Laplacian matrix. The 
application of the formula is demonstrated on three examples of networks (circuits): 
two resistive bridges and an operational amplifier. The calculated results are compared 
to results from a circuit simulator. 

In electric circuit theory, resistor networks are a special case of impedance networks. 
In general, these are defined by a complex Laplacian matrix L. Our numerical results 
indicate, that Theorem formulated in this paper holds also for systems with complex 
Laplacian matrix. Therefore, ^ should be applicable to computation of the complex 
impedance as well. 
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